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Abstract 

We consider the Navier-Stokes equation in a domain with irregular boundaries. The 
irregularity is modeled by a spatially homogeneous random process, with typical size e ^ 
1 . In the parent paper [8] , we derived a homogenized boundary condition of Navier type 
as e — > 0. We show here that for a large class of boundaries, this Navier condition provides 
a 0(e'^/^| lne|^/^) approximation in L^, instead of 0(£'^/^) for periodic irregularities. Our 
result relies on the study of an auxiliary boundary layer system. Decay properties of this 
boundary layer are deduced from a central limit theorem for dependent variables. 

Keywords: Wall laws, rough boundaries, stochastic homogenization, decay of correlations 

1 Introduction 

The concern of this paper is the effect of a rough boundary on a viscous fluid. In most 
situations of physical relevance, such effect can not be described in detail: either the precise 
shape of the roughness is unknown, or its spatial variations are too small for computational 
grids. Therefore, one may only hope to account for the averaged effect of the irregularities. 
This is the purpose of wall laws: the irregular boundary is replaced by an artificial smoothed 
one, and an artificial boundary condition (a wall law) is prescribed there, that should reflect 
the mean impact of the roughness. 

This paper is a mathematical study of wall laws, in the following simple setting: we 
consider a two-dimensional rough channel 

where = R x (0, 1) is the smooth part, is the rough part, and E = M x {0} their interface. 
We assume that the rough part has typical size e, that is 



for a ii'— Lipschitz function u; : M i— > (—1,0), X > 0. More will be assumed on the boundary 
function uj hereafter (see figure for an example of such a rough domain). 
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Figure 1: The rough domain O^. 



We assume that in this channel domain, the viscous fluid obeys to the stationary incom- 
pressible Navier-Stokes equations: 



- Au + u-Vu + Vp = 0, X e^l", 
div u = 0, X G Q^, 



Ul = <j), 

I u\dn = 0, 



(1.1) 



where denotes any vertical cross-section of and (/) > 0. The third equation in (jl.ip 
expresses that a flux (p is imposed across the channel. Note that this flux does not depend 
on the cross-section, due to the incompressibility and no-slip condition at the boundary. We 
also stress that, up to minor changes, we could apply our analysis to many variants of this 
problem, notably to elliptic type systems or to unstationary Navier-Stokes. 

In this simple setting, the search for wall laws resumes to the following problem: to find 
a boundary operator B^{x, D^), regular in e, acting at the artificial boundary S, such that 
the solution of 

- Au + u ■ Vu + Vp = 0, X 
div u = 0, X £ ^l, 

(1.2) 



[ B'(x,D,)u\ 



2 = 1 







0, 



approximates well the solution of p.ip in Q. 

This type of homogenization problems has been considered in many mathematical works. 
On wall laws for scalar elliptic equations, we refer to [2\. On wall laws for fluid flows, see 
[D El m UHl [20l [11]. See also [21| on porous boundaries. These works go along with more 
formal computations, grounded by empirical arguments (c/for instance [HI [23]). We finally 
mention [TO] for the study of roughness- induced effects on geophysical systems. 

All these studies have been carried under two assumptions: 
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• compact domains, for instance bounded channels or periodic in the variable xi. 

• periodic irregularities, meaning that the boundary function lo is periodic. 

The first restriction is just a small mathematical convenience, that gives direct compactness 
properties through Rellich type theorems. The second assumption is of course a big simpli- 
fication, both from the point of view of mathematics and physics. These assumptions were 
considerably relaxed in the recent article [8] by A. Basson and the author. As the present 
note extends this article, we now describe shortly its main results and underlying difficulties. 

In all papers on wall laws, the starting point is a formal expansion of u^: 

u^{x) ~ u^{x) + 6(t>ev{x/e) + . . . 
Formally, the leading term satisfies (|1.2p with the simple no-slip condition 

B%x,Dx)u := u = at S (1.3) 
The solution of this approximate system is the famous Poiseuille flow : 

U°{x) = {U{X2),0), U{X2) = 6(^X2(1 -X2) 

Note that is defined in all M^. This zeroth order asymptotics can be mathematically 
justified, at least for small fluxes (p: we prove in article [8] 

Theorem 1 There exists (po, Eq > 0, such that for all (j) < (po, e < £q, system (jl.ip has a 
unique solution u"^ in H^^^JQ!^). Moreover, 

\W-u^\\h\ (ne)<C^/e, \\u^-u°\\l2 (m < C'e. 

We stress that these estimates hold without further assumption on the boundary: we only 
assume that uj has values in (—1,0) and is /C— Lip. A look at the proof shows that the 
constants C and C are only increasing functions of K. 

Theorem [1] expresses that the wall law (jl.3p provides a 0(e) approximation of in 
-^u/oc(^)- ^^^o |19] for a similar result in a bounded channel. However, this wall law does 
not account for the behaviour of if near the boundary, and can therefore be reflned. Indeed, 
as the Poiseuille flow i^ does not vanish at the lower part of 50^, a boundary layer corrector 
£(l)v{x/e) must be added to the expansion. The (normalized) boundary layer v = v{y) is 
defined on the rescaled infinite domain 

1^''' = {y, y2 > c^(yi)} 

and formally satisfies the following Stokes problem 

{-A7;-FVg = 0, x&Q!'\ 
div v = 0, X G (1.4) 
v{yi:^^{yi)) = -(w(2/i),0). 

Note the inhomogeneous Dirichlet condition, that cancels the trace of . 

Although linear, the boundary layer system p.4p is quite challenging. First, the well- 
posedness is not clear. As the boundary function lo is not decreasing at infinity, one can 
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expect only local integrability of the solution v in variable yi. The derivation of local bounds 
is not obvious: the Stokes operator being vectorial, one can not use scalar tools such as the 
maximum principle or Harnack inequality. Moreover, as fi^' is unbounded in all directions, 
the Poincare inequality (which allows to get H^^^^ estimates in the channel) is not available. 
Besides the well-posedness issue, the qualitative properties of v seem also out of reach without 
further hypothesis. 

Under an assumption of periodic irregularities, the analysis of ()1.4|) becomes straightfor- 
ward. If uj is say L periodic in yi, it is easy to show well-posedness in the space 

<v£ Hl^^iflhi), V L- periodic in yi, j f \Vv\'^dy2dyi < +00 i . 
[ Jo Jujiyi) ] 

Moreover, a simple Fourier transform in yi shows that 

\\v{y)-v°°\\<Ce-^y^'^, v°° = {a,Q), a = ]- [ vi{s)ds, S > 0, 

^ Jo 

that is exponential convergence to a constant field v°° = {a, 0) at infinity. 

The constant a at infinity is then responsible for a 0{e) tangential slip. Namely, chosing 
as a wall law the Navier-slip condition 

B^{x,Dx)v = {vi — ead2Vi, V2) = at S, (1.5) 

it can be shown (in this periodic framework) that the solution of ()1.2|) provides a 0(e'^/^) 
approximation of in L^. We refer to [19j for all necessary details. The error estimate e^/^ 
comes from the fact that the boundary layer term satisfies \\e{v{x/e) — {a,0))\\i2 = 0(e^/^). 

The periodicity hypothesis is a stringent one, and one may wonder if the use of Navier 
slip condition can be justified in more general configurations. This issue has been adressed 
rigorously in the recent article [8|. Inspired by the probabilistic modeling of heterogeneous 
media (see for instance |22j). we considered irregularities that are not distributed periodically, 
but randomly, following a stationary stochastic process. Namely, the rough boundary is seen 
as a realization of a stationary spatial process. Following the well-known construction of 
Kolmogorov, this amounts to consider the space 

P = {uj-.m^ (-1, 0), ujK - Lip} 

of all possible rough boundaries, together with the cylindrical a— field C (that is generated 
by the coordinates to ^{t)) and with a stationary measure vr. Stationary means that tt is 
invariant by the group of translation 

Th P ^ P, UJ ^ U!{- + h). 

As a consequence of this modeling, the domains 0^, O''', as well as the velocity fields 
or V depend on the parameter oj. As discussed earlier, the existence result and estimates 
of theorem [T] are uniform on P. Moreover, it was shown in article [8] that the function 
to ^ u'^{uj,-) (extended by outside Q^{uj)) is measurable as a function from P to Hl^^{M?). 

Using this probabilistic structure, we have been able to extend partially the results of the 
periodic case. Key elements of our analysis are: 
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• the well-posedness of the boundary layer system, obtained in functional spaces encoding 
the relation 

v{Th{uj),yi,y2) = v{LV,yi + h,y2). 

• the convergence of v{uj,y) to (a(u;),0) as y2 +cxd, both in L'^{P) and almost surely, 
locally uniformly in yi. Such convergence is deduced from the ergodic theorem. 

More on the boundary layer system will be provided in the next sections. As regards the 
Navier wall law (jl.Sp . the main result of [8j resumes to 

Theorem 2 There exists a = a{u>) € L'^{P) such that the solution of (II. 2p . (jl.Sp satisfies 

-u^Wl^, (PxQ) = o{£). 

We remind that 1 1 If 1 1 (PxH) '•— ■^'"■^^i^ ( /p /^(a' i)nf2 l^l^^^'^^ ) 

Theorem [2] shows that a slip condition of Navier type improves the approximation of . 
As in the periodic case, the random variable a in (II. 5p comes from the convergence of the 
boundary layer v. If the measure tt is ergodic, a does not depend on a;, as pointed out in [8]. 

A natural concern about this result is the o(e) bound, which is only a slight improvement 
of the 0(e) in theorem [TJ A look at article ^ shows that this poor bound is due to the 
lack of information on the way v converges at infinity. Contrary to the periodic case, where 
convergence at exponential rate is established, the simple use of the ergodic theorem does not 
yield any speed rate. 

The present paper aims at clarifying this point. Losely, we will show that for a large 
class of boundaries, the Navier wall law provides a 0(e^/^| ln(e)|^/^) approximation of the real 
solution. Namely, we will make the two following assumptions on our random roughness: 
(HI) The measure tt is supported by 

Pa = {uJ:R^ (-1, 0), \\uj\\c2,c < Ka} 

for some a > and some > 0. 

(H2) The randon boundary has no correlation at large distances, that is the a-fields 

(T (s I— > a;(s), s < a) and a (s i— > uj{s), s > b) 

are independent for b — a > k, for some k > 0. 

Under these assumptions, the main theorem of the paper reads: 

Theorem 3 For small enough (p and under (H1)-(H2), the following refined estimate holds: 

Ik^-^^k^ (PxC)=0(e=^/2|ln(e)|V2). 

Before entering the proof of this theorem, let us give a few hints. Theorem [3] is deduced from 
a central limit theorem for the quantity v{uj,y) — (a,0). Broadly, this theorem comes from 
good properties of the random variables 

rn+l 

X^iio) = / viu;,yi,0)dyi. 

J n 
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Due to the elliptic nature of the Stokes operator, such random variables are not independent. 
However, under assumption (H2), we are able to prove that the correlation terms £'(X„Xo) 
decay fast enough as n — > oo. As a result, one can prove a central limit theorem on and 
then a similar one onv— {a, 0). We point out that such type of results for dependent variables 
with strong decay of correlations is quite classical and has been used in various fields. We 
refer to [7J for a review paper related to dynamical systems, and to recent articles \24:\ I12j for 
applications in a PDE context. 

As a consequence of this central limit theorem, we show that the boundary layer converges 

— 1/2 1 

to a constant as |?/2 I- Note that this is in sharp contrast with the periodic case, where 
exponential convergence holds (we stress that periodic boundaries are highly correlated, thus 
far from satisfying (H2)).This speed of convergence is resposible for the e'^^'^l ln(e)|^/'^ in the 
Navier wall law. 

The main difficulty is to obtain the decay of correlations of variables like X„. The proof 
relies on precise estimates of the Green function for the Stokes operator above a non flat 
boundary. Such estimates follow from sharp elliptic regularity results, where one must pay 
attention to the oscillation of the boundary. This is achieved under the regularity assumption 
(HI), using ideas of Avellaneda and Lin for homogenization of elliptic systems ^i6j. 

2 Boundary layer decay and Navier approximation 

In this section, we explain how Theorem [2] follows from estimates on the solution v of (jl.4p . 
Such estimates will be established in the following sections. At first, we remind the main 
features of v, as stated in article [8]. 

2.1 The boundary layer system 

As emphasized in the introduction, to solve ()1.4p in a deterministic way, that is for each 
possible boundary to, is still unclear. Hence, one must take advantage of the probabilistic 
setting. First, notice that a reasonable solution v should satisfy: 



Together with the stationarity assumption, this relation sort of substitutes to the identity 



used in the treatment of L— periodic roughness. It allows to extend the well-posedness result, 
through an appropriate variational formulation. 

This formulation has been described in article [8j. First, one introduces the new unknown 



v{uj,yi + h,y2). 



(2.1) 



v{yi + L,y2) 



v{yi,y2) 



w{y) := v{y) + {y2,Q)l{y^<o}{y) 



and replace system p.4p by 



< 




(2.2) 
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where [•]|y2=o denotes the jump at y2 = 0. Then, one multiphes formally the Stokes equation 
by a test function w' = w'{uj,y) that satisfies div w' = 0, w'Iq^u = 0. Integrating by parts 
over n''^n{\yi<l} yields 



/ 

Jn 



Vw ■ Vw' = w[ + {dnW — qn) w'. 

^llwil<i,w2=o) Jf7"nil?/il=n 



n>>in{\yi<i} ■J{\yi\<i,y2=o} Jn>''n{\yi\=i} 

Finally, if w, w' satisfy relation ()2.ip . one can integrate with respect to to, and thanks to the 
stationarity of vr, get rid of the annoying boundary term at the r.h.s: 



E / Vw Vw' = E 

!^"n{|j/l<l} J{|yi<l,j;2=0} 



w' 



Afterwards, this formal variational formulation can be rigorously defined and solved: in 
short, one can apply the Riesz theorem in a functional space of Sobolev type, made of functions 
w such that 

E / iVwP < +00, 



/nwn{|j(i<i} 

and satisfying almost surely ()2.ip . together with div w = 0, wIqm =0. We refer to [8] for all 
details. Note that stationarity implies: 



1 

R 



supt,RE— I \Vw\^ = E / \Vw\^ < +00. 

f7"n{|s/i-i|<R} Jn"n{|yi<i} 



Back to the original system (II. 4p . this variational solution w provides almost surely a 
solution v{uj, •) G H^oc (^''') the sense of distributions. Moreover, the ergodic theorem 
yields (see iSj) 
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15 / I^^P < +00, almost surely. 

R ^ Jn''in{\yi\<R} 

In order to understand the origin of the Navier approximation, the next step is to describe 
the behavior of v as y2 +oo. For periodic roughness, one can show exponential convergence 
of f to a constant vector field v°° = (a,0). However, the rate of convergence goes to zero 
with the period L. When dealing with stationary random boundaries, that broadly speaking 
contain all periods, the exponential decay does not hold a priori. In other words, there is a 
problem associated to the Fourier spectrum, that is discrete in the periodic case, and may 
accumulate to zero in the random case. 

Again, this problem has been (partially) overcome in [8]. The first step is to obtain a 
representation of v in terms of a Stokes double layer potential, cf proposition ??. Almost 
surely, for any 

viuJ,y) = / G{t,y2)v{uj,yi-t,0)dt 

rt 



1 r 

tdtG{y2)- I v{io,yi - s,0)dsdt 
t Jo 



where G is the Poisson type kernel for the Stokes operator over a half space. Then, the ergodic 
theorem and a few calculations yield: 



t 

f (w, yi — s, 0) ds — > w°°(a;) = (a(a;),0), t — > ±00, 
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where the convergence holds ahnost surely (locally uniformly in yi), as well as in L^{P) 
(uniformly in yi). In the case where the stationary measure vr is ergodic, the constant a does 
not depend on oj. Finally, back to the integral representation, and with similar treatment for 
derivatives of v. 

V/3GN2, |/3| > 1, E\v{;0,,y2)-a{-)\^ + yf^E 
We refer to Proposition 13] for all details. 



d^v{-,0,y2, 



V2- 



(2.3) 



2.2 Refined estimate for Navier wall law 

Most of the analysis of the present paper will be devoted to a refined asymptotic estimate of 
the boundary layer: 

Theorem 4 Under assumptions (HI), (H2), for all j3 G N^, 



(t;(-,0,y2) -«(•)) 



(2.4) 



It is of course a much sharper convergence result than ()2.3p . Before tackling its proof, we 
explain how it implies theorem [21 Arguments are direct adaptation from section 5 in [8j. 

On the basis of the boundary layer analysis, one can build an approximation of of 
boundary layer type. Namely, we introduce 

u^pp(w,x) = vP{x) + ^(pev {uj,—^ + Q4)ev}'{ijj,x) + Q(l}er^{uj,x). 

In this approximation, is the Poiseuille flow and v{ijJ,-) is the boundary layer solution of 
(|1.4p . As V does not converge to zero at infinity, we add a large scale corrector v} satisfying: 



^ ■ Vu^ + ■ V-u^ - Au^ + Vp = 0, X £n, 
div = 0, x G 0, 



/ ■ eidx2 = 0, 
Jo 



(2.5) 



U^\y2=0 = 0, 



y2=i 



-(a,0). 



It is just a combination of a Couette and a Poiseuille flow: v} = ax2 (2 — 3x2) ei. Still, this 
approximation does not vanish at the boundary, which explains the addition of another term 
r^(u;,x). It must satisfy 

r^(xi,0) = 0, 



r^(xi,l) = vl"^,-] -(a,0). 



£ e 

^ div = 0, X G $7. 
This remainder can be taken small in the sense of 

Proposition 5 This problem possesses a (non unique) solution r^ such that 



(2.6) 
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Proof: The proof of this result mimics the one of proposition 14 in [8]. The corrector can 
be chosen in the form 

= V^V) V' = cl{xi)x2 + b{xi)x2 + c{xi)x2 + d{xi). 

The streamfmiction ip is determined up to a constant and polynomial in X2- Its coefficients 
have explicit dependence on v — (a, 0). Hence, the estimate on r'^ follows from the control 
of various terms involving v — (a, 0). For instance, one must bound the L^{P x (—1, 1)) norm 
of 



^1 / t 1' 
V2 { uj,-,- ] dt 



V2 iuj,yi, - ] dyi 

l/e 



{vi - a) [uj, — ,^2] dy2 

lui(xi/e) V £ / 

'■l/e 

{vi - a){uj,0,y2)dy2 :=/^(a;,3;i) - I^{uj,0) 

where the last equality comes from the Stokes formula. Using stationarity of vr, we get 

f^i / t l\ ^ 
E xi^ V2(-,-,-)dt < 4E|r(-,0)p. 

Thanks to the refined estimate (j2.4p . we finally obtain 



E|/^(-,l 



< C [e^E 



oj{0) 



l/e 



\{vi-a){-,0,y2rdy2 +e / E\{vi - a){-,0,y2)\'dy.. 



rl/e 

< C'e'^ + C"e j y^^ dy2 = 0{e\ lne|) 



All other terms involve similar computations. These are straightforwardly adapted from the 
proof of proposition 14 in [8j, using (j2.4p instead of (|2.3p . 

Once the approximate solution n^^p is built, one can obtain by energy estimates the 
following bounds, for cp small enough: 



Af II 



0(e3/2|in(e)|V2 
: 0(e3/2|in(e)|V2 



which of course imply theorem [2l As the proof is very similar to what was done in paper [8], 
we do not expand more and refer to it for all details. 



3 A central limit theorem 

Up to the end of the paper, we will assume (H1)-(H2), and focus on theoremfH It is classical 
that (H2) implies ergodicity of tt, so that the constant a does not depend on uj. We start 
again from an integral representation 

a^(t;(a;,0,y2)-(a,0)) = \ d^' d^^^ G{t,y2) {v{uj, -t,0) - {a,0)) dt 

t (3-1) 

= - I d^'-^'d^^Git,y2) [ {v{u;,-s,0)-{a,0)))dsdt, 

JR Jo 
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where the matrix kernel G is given by 

We introduce ^ 

V{io,t) := / {v{u,-s,0) - {a,0)) ds. 
Jo 

A simple change of variable leads to 

yf+^/\^7;(^,0,y2)-(a,0)) = / dt'^'d^'^G{u,l) y~'/' V{co,y2u) du. (3.2) 

JR 

Our first goal is to show that the l.h.s. converges in law to a gaussian distribution, for all /3. 
We will focus on the case = 0, the other cases being handled in the exact same way. We 
state 

Proposition 6 The function V satisfies the following properties: 

i) E|y(-,t)|2 < c\t\ 

ii) The random process ^("^i ^2 u) converges weakly to a gaussian process B{uj,u) as 1/2 

goes to infinity. 

iii) The covariance matrices also converge, that is for all indices i,j and for all s,t, 

E 2/2^1 V,{;y2 s) Vj{-,y2 t) > EBi(-, 8)B,{-, t) 

J/2-++00 

We remind that the process X'"'{ijj,t) with values in converges weakly to X{uj,t) if, for all 
T > and all continuous bounded function : C ([— T, T],M^) 1— > R, 

EJ^(X") > EJ^(X). 

n^+00 

Theorem is then a direct consequence of 

1/2 

Corollary 1 The random process 7/2 (^('-^) 0,7/2) — («)0)) converges in law to a gaussian 
vector with zero average. Moreover, for all i,j, 

EK(-,0,y2) - (a,0)i) {vj{-,0,y2) - {a,0)j) > aij, 

where a is the covariance matrix of this gaussian vector. 

Proof of the corollary: To prove convergence in law to a gaussian vector Ma of covariance 
matrix a, we need to show that for any F G C^(M), 

Ef( [ dtG{t,l)y-^/^V{-,y2t)dt) >EF{Ma). 

Unsurprisingly, we take 

K ■■= [ dtG{t,l)B{uj,t)dt. 
Jr 
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We decompose, for any T > 0, 



E F (^J^ dtG{t, 1) V{:,y2t) d?j-EF (AA,) 

y ^ dtGit, 1) ^2^'/' 2/2t) dtj-EFU ^ dtGit, 1) t) dtj 
+ E F atG(t, 1) ^2^'/' yst) dt) - E F (^^^ atG(t, 1) ^2^'/' ^2*) d?j 
+ EF^y" dtG{t,l)B{-,t)dt^ -EF (^j dtGit, 1) B{-,t) dt 



:= Ji + J2 + J3 

We now show that expressions Ji, J2, converge to zero (J3 is similar to J2 and simpler). Let 
6 > 0. We have 



1/2 / „ \ 1/2 



IJ2I < max|F'|E [ \dtG{t,l)\ y^^^^ V{uj,y2t) 

J\t\>T 

< max\F'\ ( [ \dtG{t,l)\dt] ([ \dtG{t,l)\ E [y^^ \V{uj,y2t)\^ ) dt 
\J\t\>T J y\t\>T / 

<C\\ \dtG[t,\)\dt\ If \dtGit,l)\tdt] 

\J\t\>T J \J\i\>T J 

where the last line comes from point ii) of proposition [6l Thus, for T large enough, indepen- 
dently of y2, I J2I < '^/2. Such T being fixed, for y2 large enough, we get | Ji| < 5/2 by point 
i) of proposition El This yields convergence in law. The convergence of the covariance matrix 

E{vi{;0,y2) - (a,0),) {vj{;0,y2) - (a,0),) 

= / [ E(dtGis,l)y^'/^Vi;y2s)) (dtG{t,l)y-'^^ Vi;y2t)) dsdt 

follows from the dominated convergence theorem, using i) and iii) of proposition El We get 

a^j = E [ [ {dtG{s,l)B{;s)), {dtG{t,l)B{;t)) . dsdt. 
Jr Jr 

This concludes the proof of the corollary. 

It remains to prove theorem El Note that point ii) is essentially a central limit theorem 
for the sequence of random variables 



X'^iuj) = F oTniuj), F{uj)= ! {v{uj,t,0) - {a{uj),0)) dt. 

Jo 

The problem is that these random variables are not independent, due to "propagation of 
information at infinite speed" in the Stokes system. To establish a central limit theorem for 
such type of sequences is a classical question. The basic idea is that one can extend the central 
limit theorem to non independent sequences that feature a good decay of correlations as n goes 
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to infinity. We now illustrate this general principle on our problem, using assumption (H2). 
We follow the presentation of article [21], in which a similar question arises for a semilinear 
heat equation with random source. Let the cr-algebra generated by the applications yi i— > 
^{yi)j \yi\ < n.We state the following lemma: 

Lemma 7 Suppose that := E(f(-,0,0) | C") satisfies 

Elv" -t;(-, 0,0)1^ < Cn-" 

for some a > 1. Then, proposition\^ holds. 

Proof of the lemma: We write the decomposition 

+ 00 



f(-,0,0) - (a,0) = - (a,0) + 



2J ^2^-1 



with the sum converging in L?'{P). The corresponding sum for V is 



+ 00 

V = V V^{u;,t)= / 

,=0 -^0 



2J 



V I o Ts{uj) ds, 



where v^/^ (q,^q)_ xhen, we have: t)||i2(p) < Xlj!^ t)||i2(p). By the as 

sumption of independence at large distances, the correlations E oTt{uj)v'^^ o Tt+s{'-^ 
and E (v'^^ o Tt{uj) v'^^ ^ o r(-(_s(cj) ) vanish for \s\ > k + 2^'^^. We introduce 



n :- 



[|t|/(^ + 2^-+i)]. 



If re = 0, we just write 

If re > 1, we decompose 
E\V^{-,t)\ 



E\VH-,t)\ < ItrE 



2J 2^~^ 
W — V 



E 



< E 



r(k+l)t/n 



fc=0 



kt/n 



f o — V oTg] ds 



j-t/n 

^0 S 



°Ts+kt/n-V^^ °Ts+kt/n 



ds 



< 2(^ + 2^+1) / Ve 
■^0 fc=0 



23 2^-1 
f — f 



Using the bound on the conditional expectations, we end up with 

\\V\;t)\\hiP) < C\t\ min(|t|,2^-)2-^°. 



(3.3) 



for some constant C = C{k). We thus get i). To prove ii), we just write the decomposition 



;(-,0,0) -(a,0) = 

3=0 



23 2^-1 
V — V 



, y-^'^V{u:,ty2) = ^'(^'*2^2). 

3=0 



-1/2 ' 



+ 00 
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It is well-known that each finite sum satisfies a central limit theorem, that is 

k 
3=0 

in the sense of weak convergence, to some gaussian process B^{uj, t). Moreover, the covariance 
matrix also converges, that is 

k k 

y2-iE^ V,\;sy2) J] VH^ty^) — - ESf (u;, s) t). 

3=0 3=0 

In short, it is due to the fact that the random variables 

X'^'^uj) = o r„(cj), F\uj) = (^v"^' - v'^'~^^ o Tt{uj) dt, n e Z, 

have zero correlations at large distances: see |13i theorem (7.11) p424] for a similar result and 
detailed proof. Moreover, thanks to estimate p.3p . the remainder 

+ 00 

j=k 

converges to zero as /c — > +cxd, locally uniformly in t, uniformly in y2- Hence, points (ii) and 
(iii) of proposition [6] hold, which ends the proof of the lemma. 

We still have to estimate the variance of — v{-,0,0). Following [23], we can turn this 
question into a question of domain of dependence for solutions of (jl.4p . Precisely, starting 
from the measure vr on P, we define a new measure tt"' on the product space 

= {{0Jl,0J2) ePx P, UJl{t) = UJ2{t), \t\ < n} . 

endowed with its cylindrical cr— field. Namely, vr" is defined in the following way: 

1. 7r"(yl X A) := 7r{A), VA e C", which determines vr" over the sub cj— field generated 
by the applications 1 1-^ {uji{t),uj2{t)), \t\ < n. 

2. For ah k > 1, for all t^,...,t'' with \P\ > n, for all borelian subsets Bl,...,Bi, 
B^,...,B^ of M, and 

Ai := n';^-^{uji, uJi{tj) eBi), A2 := n^^=i{a;2, Mh) ^ 

7r"(yli X A2 I W2) := vr(Ai | 7r(^2 | C")(u;2), which determines tt" condi- 

tionaly to P". 

It is then easy to derive the following identity, see [24j : 

E|i;"-T;(-,0,0)|2 = i / 0, 0) - ?;(w2, 0, 0)|' 

Thus, if and Q''\uj2) are two boundary layer domains with boundaries that coincide 

over [— n, n], we need to estimate the difference of the corresponding boundary layer solutions 
u(a;i,0, 0) and f (012, 0,0). This is the purpose of the next section. 
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4 Decay of correlations 



Throughout the rest of the paper, we will assume (HI). The main difficuly is to prove the 
following 

Proposition 8 Under assumption (HI), for all < t < 1, for almost every uji,uJ2 S P, 

\v{u;i,0, 0) - v{oJ2, 0, 0)1 < (4.1) 
where C does not depend on uji,uj2- 

Together with the results of the preceding section, this proposition concludes the proof of 
theorem m (take r > 3/4), and therefore the proof of the main theorem[2l In fact, the sharper 
bound 

\v{ui,0,0)-v{uJ2,0,0)\ < -, 

n 

that is with r = 1 would still be true. We will discuss this briefly in the last section of the 
paper. For the sake of brevity, we only prove here the weaker form (j4.ip . 

The main difficulty is that the boundary layer solutions v{uJi,y) and v{uj2,y) of ()1.4p are 
not defined on the same domain, so that estimates on the difference are not directly available. 
If the Poisson equation rather than the Stokes system was considered, representation of the 
solution in terms of Brownian motion would allow to conclude quite easily. Again, this will 
be explained in the last section of the paper. 

In the case of system (|1.4p . we are not aware of such representation, and the bound (j4.ip 
will come from an accurate description of the (matrix) Green function of the Stokes operator 
above a humped boundary. We consider for all uj G C^'", and for all z G fi'''(a;) = {y2 > 
cj(yi)}, the system: 

( -AGUz, •) + VP^(z, •) = 5, h in 17^'(u;), 

divG^(z,-) = \n^^\uj), (4.2) 
G^(z,-) = oiidQ!'\uj). 

where 5z is the Dirac mass at z, and I2 is the 2x2 identity matrix. Let us remind how to 
build the matrix Green function {Gi^,P^^). Up to a vertical translation of the domain, we 
can first assume that Z2 > 0. We then introduce the Green function (Go,Po) for the Stokes 
operator in the upper-half plane, see [H]. Extending Go{z, ■), Po{z, ■) by for y2 < 0, the 
functions 

H{z, •) := GUz, •) - Go{z, •), Q{z, •) := P^z, •) - Po{z, •) 

satisfy formally 

' -AH{z,-) +VQ{z,-) =0 inQ^'^uj), 
dwH{z,-)=0 mn''\uj), 
H^{z,-)=0 ondn^^iuj), 

[H{z, •)] = 0, [d2H{z, •) - Q{z, •) es] = - [d2Go{z, •) - Po{z, •) 62] , 
where [ • ] is the jump along {ya = 0} n n^^uj). The jump on the derivative is explicit, as 

d2Goiz, (,„ 0-)) - Poiz, (,„ 0-)) ^ e. = (/f -f - f )^^) . 

tt{{zi - yiy + z^y \[zi - yi)y2 ^3 / 
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Standard variational formulation yields existence and uniqueness of a solution H(z, •) with 
VH(z, •) in L^. In turn, this provides a unique solution Gi^{z, •) to (|4.2|) . The corresponding 
pressure field Pi^{z,-) is determined up to the addition of a constant matrix. Note that 
uniqueness yields the relation 

Gr^X'^){z,y) = G^{{zi + h,Z2),{yi + h,y2)). (4.3) 

Our key estimate is provided by 
Lemma 9 For all < t < 1, for all z,y ^ satisfying \z — y\ > \, we have 

Y: \d^yGUz,y)\ + \VyPUz,y)\ < c '^(^ra + y)\ (4.4) 

where 6{-) denotes the distance to the boundary dVl!'\uj), and C is a constant depending only 
on T and on \\oj\\Q2,a. 

Note that by symmetry of G, we also have 

|/3|<2 ' ^' 

Moreover, in the course of the proof of lemma [H we will show that for all G VI^^lx)), 

\G^{z,y)\ < C (|ln|z-y|| + l). (4.6) 

Let us show proposition El postponing the proof of lemma to the next section. We first 
need to connect the solution v{uj, ■) of p.4|) to the Green function G^^. For this purpose, we 
rather consider 

w{u},y) = u(u;,y) +2/2l{J;2<o}(y)■ 
which satisfies ()2.2p . Note that v and w coincide at y = 0. Formally, w should be equal to 



wyijj, z^ 

'{3/2=0} 



/ G^{z,y)eidy. (4.7) 



Using estimates (|4.4|) . (|4.6p . it is standard to show that w is a, solution of in H}^^{Q.u). 
Using bound (14. Sp . one has even 



iVwP < C < +CX), 



/n"n{|2i|<i} 
for all u; G Pq. 

Extending w{u}, ■) by outside ^'''(u;), one can show that u) ^ w{uj, •) is measurable from 
Pa to HI^^{M?) (see the appendix for details). Moreover, thanks to (14. 3p . w satisfies the 
stationarity relation w{Th{uj),y) = wiuj, {yi + h, 2/2))- Finally, using that w and w both satisfy 
(12. 2p . a simple energy estimate on the difference leads to 

E / \V{w-w)f = 

Jn>"n{\zi\<i} 
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which shows that w = w almost surely. 
It remains to estimate the difference 

t;(a;i,0,0) -v{uj2,0,0) -- 



(G^,((0,0),y) -G^,((0,0),y)) ei, 

{2/2=0} 

for every coi, uj2 in Pa which coincide over [— n, n]. This integral is bounded by 



h+h 



J/2=0,|?/i|>n 



\G^,-G^,\iiO,0),y)dy 



+ 



J/2=0,|j/i|<n 



\G^,-G^,\ {(0,0), y)dy 



The use of (|4.4p gives 



h < G 



y2=0,\yi\>n \yi 



1 , c 



where C, which depends a priori on ||a;j||(72,Q , can be chosen uniformly over Pq^ 5 ^h G ' 
norms are bounded by K^- To bound the second term, we first assume that uj2 > for 
\yi\ > n, which is always possible up to introduce an intermediate third boundary. Hence, 
Q,^''{u)2) C Q^^LJi). To lighten notations, we introduce 



as well as 



P{y) := P^,((0,0),(yi,cu2(yi))), V 



which defines a continuous extension of Puj2{{0,0), ■) outside Q^''{ll!2)- Finally, we define the 
vector fields 

U{y) := {G^,-G^,) {{0,0), y), Q{y) := {P^, - P^,) {{0,0), y), y G J^^'(t^2), 
U{y) := G^,{{0,0),y), Q{y) := P^,{{0,0),y) - P{y), y e n\\2- 

They satisfy 

' -AU + VQ = 0, ye 9^\u)2), 
-AU + VQ = -VP, yG05^2, 
divC/ = 0, yGO^'(wi), 
U = 0, yeda^\uji), 
[^]|ri,2 = 0, [a„f/-Q®n]|r,2 

A direct energy estimate yields 



(4i 



|5„G^2((0,0),y)||C/| + 



1,2 



''l,2 



-5„G^,((0,0),y)|r,,2 
|VP| \U\ 



< 



1/2 



< G 



a„G^,((0,0),y)p 



1/2 



+ 



ivpr 



^'1,2 



1/2 



''l,2 



1/2 
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Note that all y in both ri_2 and satisfy \yi\ > n. Using (j4.4p . we end up with 

/ \VUf < Cn^-^^ 

Back to I2, we obtain 

This ends the proof of proposition ([8]). 

5 Green function estimates 

This section is devoted to the proof of lemma\^ that is sharp estimates on the Green function 
for the Stokes operator above the humped boundary 2/2 = "^(yi); where uj belongs 
to C^'". A fundamental remark is that the Green function satisfies the scaling 

Ve > 0, G^e{ez,ey) = G^{z,y), = ea;(xi/e). (5.1) 

We want estimates ()4.4p to hold for — y| large, that is for e := \z — y\~'^ small. By 
relation (j5.ip , to establish such estimates amounts to get local estimates for the Green function 
Gf^e. Thus, this is again a homogenization problem: more precisely, we must show that the 
oscillations of the boundary at frequency do not affect too much the estimates on G^^e, 
so that it behaves as the Green function for a half-plane. A very close problem has been 
considered in the papers [5l[6] by Avellaneda and Lin, namely the derivation of local estimates 
for elliptic systems div (A(x/e)V-), in which j4 is a positive definite matrix with periodic 
coefficients. Our reasoning follows these papers. 

For all X G M^, r > 0, we will denote -D(x, r) the disk of center x and radius r, and 

D^{x,r) := D{x,r)n{x2> uj^ixi)}, T^{x,r) := D{x,r) D {x2 = uj^{xi)}. 

An important property is that for all < r < 1, 

\D'{x,r)\ > r]r'^, (5.2) 

for some r] > independent of e. More precisely, r] only involves the Lipschitz norm of uj^, 
wich is bounded uniformly in e. This will be used implicitly throughout the sequel. 

The core of the proof is to derive elliptic estimates uniform with respect to e on the 
following Stokes problem: 

' -An + Vp = div /, x £ D%xo,l) 

divn = 0, xeD^xoA), (5.3) 
n = 0, x£r'{xo,l), 
where xq € M^. More precisely, there are two steps in the proof: 

1. We show a e-uniform Holder estimate on u: for all / G L*?, q > 2 and for /_i = 1 — 2/q, 

\W\\cO:t^{D^{xo,l/2)) < C {\\f\\Li{D-'{xo,l)) + \W\\L^{D^ixo,l))) ■ (5-4) 
where C depends only on I I I I . 

2. Thanks to this Holder estimate, we prove (|4.4p . 
The two next paragraphs correspond to these steps. 
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5.1 Holder estimate 

To obtain a Holder regularity result, a classical approach is to use a characterization of Holder 
spaces due to Campanato (see [Ej): for 17 an open connected bounded set, v E C°''^(0) iff 
V £ L^(r2) and 

sup ]^ I \v -Vx,r\'^ < CO, il{x,r) := Qn D{x,r), Vx,r ■= ,„/ / v. 

xGn,r>0 Jn(x,r) M^, r)\ jQ(^x,r) 

One then tries to control such local integrals through energy estimates. This approach has 
been successful in the study of elliptic systems, see the work of Giaquinta and coauthors [E]. 
It extends to the Stokes type equations, c/ article [T7]. For us, it amounts to controlling 

■^x,r ■~ 94_9,, / 1^ ~ r\ < OO, Ux r '■— 77^77 \T / ^ 

where u is solution of (|5.3p . Note that, thanks to (|5.2p (see [16]). 

\\u\\co^^^{D^ixo,i/2)) < Cxo \\\u\\l2{d^ (xo. 1/2)) + sup I%x,r)\ 

\ xeD''{xo,l/2),r>0 J 

with CxQ independent of e. In our case, the main problem is to keep track of the dependence 
of 7| on e. It involves a discussion of the way e relates to r. Broadly speaking, the idea is 
the following: if r is large compared to e, then the oscillations have small enough amplitude 
to apply the regularity results of the flat case. On the contrary, when r gets as small or 
even smaller than e, one can rescale everything by a factor e, so that oscillations of the 
boundary have frequency 0(1), and are no longer annoying. Implementation of this idea is a 
bit technical, and follows closely the work of Avellaneda and Lin. 

We first remind a few elements of regularity theory for Stokes type systems. Let an 
open connected bounded set, with Lipschitz boundary. Then, for any (p G L^(J7) satisfying 
J^ip = 0, the problem 

div w = if, w\dn = 

has one solution w satisfying < C* ||99||^2(q), where C can be taken as an increasing 

function of and of the Lipschitz constant K of the boundary, see [14] . Thanks to this 
result, one has quite easily, 

\\P- I p\\l-^{q.) < C'||Au + /||||;^-i(f^) (5.5) 
Jn 

where {u,p) G H^{Q,) x L^(r2) satisfies (in the distributional sense) 

- Au + Vp = f, div u = 0, X £ n. (5.6) 

Again, the constant C in ()5.5p depends only on and the Lipschitz constant of the boundary. 

We now state the famous Cacciopoli inequality: 
Lemma 10 For all < r < 1, any solution u £ H^(Q) of (|5.3p satisfies 

W'^uWl^D^ix^r)) < C (r~^ \\u\\L^DHx,2r)) + W f W Li {D^ {x ,2r))) ■ (5-7) 
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Sketch of proof: We remind the main elements of proof. Let r/ a smooth function with compact 
support in D{x,2r), with = 1 on D(x,r). Note that |Vr/| < Cr^^. Multiplying (j5.3|) by the 
test function r/^ti, and integrating by parts, one has easily 

[ |Vn|2 < [ r/2|Vn|2 < C r'^ [ |n|2 + C [ |/|^ 

JD'^{x,r) JD'^{x,2t) JD'^{x,2r) J D'^{x,2r) 

+ \\p-Vx,2r\\L'\D^{x,2r)) l|div 11^2 (£,s(^^2r)) • 

Using ()5.5p . we get 

Hp -Px,2rllL2(D^(x,2r)) ^ C HAu + div /||/f-l(D^(x,2r)) = C|b II Hi (D^ (x,2r)) 

where v E Hq{D^{x, 2r)) is the solution of 

— Af + Vp = Au + div /, div V = 0, 'y|a_D'=(a;,2r) = 0- 

Note that the previous bound is uniform in e, as it only involves the Lipschitz constant of ui^ 
which is uniformly bounded. A simple energy estimate on v gives 

l|Vu||L2(£,£(a,^2r)) < C (l|V'"llL2{D=(a:,2r)) + ll/llL2(D'^(x,2r))) 

As div (r]u) = Vr/ • u, and using Holder inequality on /, we end up with 
f iV^xp < /" rf\Vu\'' < Cr-^ I |«|' + r^A' 

JD^{x,r) JD^{x,2r) JD^{x,2r) 

+ '^ll^'"llL2(D=(a;,2r))' 

where (5 > is arbitrary small. We conclude as in [171 Theorem 1.1, page 180]. 

Inequality of type (|5.7|) has been used by Giaquinta and Modica in the study of elliptic 
regularity. In the context of Stokes type system, they obtain a local estimate, see |17] : 



Theorem 11 Let of class , and {u,p,f) G H^{Q) x L^(il) x L'^{i^), q > 2, satisfying 

-Au + Vp = div f, div u = 0, xeO(xo,l), wbnnDCxo,!) = 0- 
Then, u G C°'^'{n{xo, 1/2)) for = I - 2/q, and 

\\u\\cO't^iQ{xo,l/2)) < (ll'"llL2(n{a;o,l)) + ll/llL9(n(a:o,l))) ■ i^-^) 

Unfortunately, we cannot use this theorem assuch. Indeed, the constant C in the last reg- 
ularity estimate involves the modulus of continuity of V7, where X2 = j{xi) describes the 
boundary. In our case 7 = 0;'^, such modulus of continuity is not uniformly bounded in e. We 
must proceed in several steps to control the local integrals Note that theorem 1111 implies 
estimate (j5.4p when D'^{x(), 1) is far from the boundary. Thus, we can restrict ourselves to a 
case in which xq is close to the oscillating boundary, for instance belongs to the axis X2 = 0. 

Lemma 12 For all 9 small enough, there exists Eq > such that for all e < Eq, and for all 
solutions of ([531) satisfying ||u^||L2(£,s(a,(,^i/4)) < 1, ||/||l'j(d^(xo,i/4)) < ^o, one has 



P llL2{D=(xo,e)) 
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Proof of the lemma: Suppose that the result does not hold. Then one can find 9 arbitrary 
small, and sequences u'^^ , f^ satisfying 

ll'"^^llL2{D'=i(xo,l/4)) - 1' ll/'' llL2(D=i(xo,l/4)) . , * 0' IW'^ W {D' J (xo ,9)) > 

One can extend all the Uj , /-^ by outside D^^ (xq ,1/4) so that all these functions are defined 
on the fixed domain D{xo, 1/4). From the bound on u^^ , up to extract a subsequence, we 
get 

n^^ converges weakly to some u in L'^(D(xo, 1/4)) 
and by Cacciopoli inequality (|5.7|) . 

n'^J converges weakly to u in H^{D{xo, 1/8)), and strongly in L?{D{xq, 1/8)). 

One can then take the limit in (15. 3p . which yields 



-Au + Vp = 0, div u = 0, in L>(3;o,l/8) n {x2 > 0}, ^i|D(xo,i/8)n{x2=o} = 0- 

As the upper half plane is a regular domain, we can apply theorem II 1|, so that for all /!>//, 
for ah e, 

lkllL2(D(xo,e)n{a;2>0}) < ^TT ||n||cO,A(l3(a;o,6»)n{x2>0}) ^''^"^ 

< ll'"llL2(D(xo,l/8)n{a;2>0}) ^'"'^^ < C O'^^^ . 

For 9 small enough, it contradicts the lower bound on \\u'^^\\L-2(^j:)<^j (^xofi))- 

We fix 9, £q as in lemma [TJJ We then state 
Lemma 13 For all £,k satisfying e/9'^ < eo (k>l), 

2 



, , < 0'^k^l+2 M|'U^||^2(£,a:(0,l/4)) + — II/IIl2(£,.(0,1/4)) ) (5. 



Proof of the lemma: The lemma is deduced from an induction argument on k. For /c = 1, the 
bound (|5.9|) is given by lemma WI\ Assume now that this bound holds for A: > 1. Up to a 
horizontal translation, we can assume that xq = 0. Then, we set 

1 

M := ||n^||L2(D^(0,l/4)) + — II/IIl2(/3e(o,i/4)), 

^0 

and introduce the rescaled functions 

V := 9-^^'M~^u{9^x), g = 9^-^^'M-^ f{9^x). 

They satisfy 

-Aw + = div 5, div u = 0, X e -D^(0,6'~''/4), ?^|r='(o,e-fc/4)- 
Moreover, one has ||/||l'?(Z)=(o,i/4)) ^ ^O) and by the induction assumption 

II'^IIl2(Z)=(0,1/4)) < 1- 

Applying lemma [T2] to v and g yields the result. 

We can now finish the proof of estimate (|5.4|) . Let x G D'^{xq,1/2). We need to bound 
P{x,r), for r > 0. There are two cases to handle: 
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The distance between x and the boundary {x2 = uj^{xi)} satisfies 6^{x) > 



eo 



Up to take a smaller £q, we can suppose that ^ > e, which implies that there exists x'q on 
the axis {x2 = 0} with | I < 6^{x). By lemma [m for all e/eo < r < 1/12, 

[ |n^|2 < Cr^^^+^hu'U^f^nHx'oA/m + \\f\\LHDH^'o,m)Y ■ (5-10) 

If r > 6'^{x)/2, D^{x,r) C D^{xQ,3r), and the previous line implies 

|^£|2 < Cr^''+^ (\\u'\\L2^D^^x'm) + \\f\\Li{DHx'm)Y- 

On the contrary, if r < 6^{x) /2, then D^(x, r) = D(x, r) (it does not intersect the boundary). 
A simple rescaling of ()5.8p yields 



\\u\\cO-i'iD{x,5''{x)/2)) < C[6'^{x) ^ '^\\u^\\l'2{D{x,S'^{x))) + II/IIl9(D(x,<5=(x)))) • 

Thus, 



JD^{x,r) 



Now, by lemma [T3| as > e/eq, 

ll^^llL2(D(a',5^(a;))) < ll""-^ IIl2(£i(^^_2<5'=(x))) 

< C5^{xY+^ (||'"^IIl2(D-(xJ„1/4)) + II/IIl9(D-(x[„1/4))) • 

Using the two last inequalities, we end up again with 

/ \u^\^ < Cr^''+^ (\\u''\\L2^D^^x'm) + II/IIl9(D-K,1/4)))\ 

JD^{x,r) V " " / 

which in turn clearly implies 

/ 1^^^ -^^X,rP < C'r2^+2 ('||n^||i2(£,.(^/^^l/4)) + 1/4))) , 

JD^{x,r) 

As D^{xq, 1/4)) C D'^{xq, 1), this gives the required estimate. 

• The distance between x and the boundary {x2 = lo^{xi)} satisfies 5^{x) < 



e 

eo ■ 



This time, there exists Xq on the axis {x2 = 0} such that jx — Xq| < (5^(x) + e < 2e/eo. Again, 
for all e/eq <r< 1/12, D{x,r) C D{xQ,3r) and (jS.lOp implies 

D'^{x,r) ^ ' 

It remains to handle the case r < e/eq. Up to a horizontal translation, we can assume that 
Xq = 0. We introduce the rescaled functions 

— ) (—^ ' 9 = f — ) / (—X 
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They satisfy in particular 

- Av + Vq = diY g, div v = 0, x e D^°{0,1), t;|r.o(o,i) = 0. (5.11) 

It is a Stokes type system set in a domain independent of the small parameter e. Hence, 
we can apply theorem II 11 which yields: for all r < 1, 



/ b-l'x.rP < C'||w||cO.M{D=o(x,r))?'^''^^ < C (11^^1^2(13=0(0,2)) + \\9\\l^{D-'0{0,2))) 

JD^O {x,r) 
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Back to the original unknowns , f, we obtain the control of /| for r < e/eq. This ends 
the proof. 

5.2 Bounds on the Green function 

From the above Holder estimate, we can deduce the estimate (14. 4p . The arguments are again 
adapted from article [5l pages 819,829-831]. For the sake of completeness, we describe the 
ideas at play. The first step is to establish the following bound: for all x,x' in {x2 > cij^(xi)}, 

\G^.{x,x')\ < C7 (|ln|x-x'|| + l) . (5.12) 

where C only involves ||a;||(7i,Q. Note that it implies ()4.6p . To lighten the notations, we 
drop the suffix lo, denoting G^, G instead of Gi^e,Gi^. Let us introduce the Green function 
G^{x,t,x' ,t') for the Stokes operator over {x2 > lo^{xi)} x T. Namely, it satisfies for all 
(x,t) G {X2 > io'ixi)} X T 

( -AG'{x,t,-)+VP'{x,t,-) = 6x,th in{x2 >cj^(xi)} xT, 

div ^^(x,*,-) = in {x2 > uj'ixi)} xT, (5.13) 
G%x,t,-)=0 on {x2 = uj%xi)} xT. 

One has easily that 



G^x,x') = j (Gf(x,0,x',t'),G|(x,0,x',t')) dt' . 
The point is to show that 

\&{x,t,x\t')\ < C- n^-T--n- 

\x — x'\ + \t — t'\ 

The estimate (|5.12|) is then obtained by integration with respect to t', at t = 0. Such bound on 
will be deduced from a repeated use of the Holder estimate (15. 4p . Note that this estimate 
extends without difficulty to similar Stokes problems in dimension n > 2, with q > n and 
/i = 1 — n/q. In particular, it holds when the domain is {x2 > uj^{xi)} x T. 

Namely, let x := {x,t) and x' := {x',t'). Let r := \x — x'\, and / G C^{D'^{x' ,r/3)). We 
consider the quantity 



u%i) = / G'{x,2)f{z)dz 
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The field satisfies a Stokes equation witli source term / over {x2 > lu^{xi)} x T, with a 
Dirichlet boundary condition. We therefore apply the estimate (|5.4|) to u^. Properly rescaled, 
it yields 

\u'{x)\ < C [ \u'\A , 

\ JD'^(x,r/3) J 



where we used the fact that / vanishes over D^{x,r/3). Thus, we get 



DE(i',r/3) 



G'ix,2)f{z)dz 



1/2 



,£|2 



< C I r " / 

D^(x,r/3) 

\ 1/6 

u'f] < C'r-^l'^ I / \U 

D^{x,r/3) J \JD^{x,r/3) 

\J{x2>ui^ixi)}xT J 



1/6 



e|6 



< C"r'/^f\\%' 



where the last two inequalities come respectively from the Sobolev imbedding theorem (note 
that is zero at the boundary so that the imbedding does not involve lower order terms), 
and from the standard energy estimate on the Stokes system. By duality, we infer that 



,-3 



\ 1/2 

\&(x,z)\^dz] < Cr-^. 



_De(£',r/3) 



Using that G^{x,-) satisfies a homogenenous Stokes system over D(x' ,r /3), one more appli- 
cation of ()5.4p leads to 

&{x,y) < C / \G^x,z)\'^dz\ < C r-\ 

\ JD^{x',r/3} J 



Inequality (|5.12|) at hand, we can derive the final estimate on G. Let x, x' £ {x2 > uj^{xi)}. 
Set this time r := \x — x'\. For all x such that |x — x| < 2r, (15.121) implies 

\G^x,x')\ < G {\ln\x-x'\\+l) < C"(|ln|rc-x'|| +1). 

Applying ()5.4p to the function we get for any r € (0, 1) 



\G%X,X')\ < d%xY\\G%;x')\\cO,.!^nix,r/3)) 

< GrS'{xYr-^-^\\G%;x' 



that leads to 



\G'{x,x')\ < Gr (|ln|a;-x'|| + 1) 



L2(D(s,2r/3)) 

6%xy 



\X — X' 



Now reversing the roles of x and x' , we obtain 

\G%x,x')\ < Gr (|ln|a;-x'|| + i; 



6%xY6^{x'y 



x — x 



I\2t 
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Using the scaling relation ()5.ip . we get for all y, z G {1/2 > Lj(yi)}, for e := \y — z\ ^, 
a(.,y) = G'(s.,ey) < (Iln - ,|) I + 1) ^"'^^ ^''^5 = c/'''^ ''^'^ 



|e(y-z)|2^ \y-z\ 



2t 



Using classical local regularity results for the Stokes equation in a C^'" domain (see ^l7\ 
Theorem 1.3, page 198], which extends theorem 1 11|) : for I2: — ?/| > 1, 

^ \d^Giz,y)\ + \VyP{z,y)\ < C ||G(z, •)||L2(Dfe,i/2)) 

l/3|<2 

5{zy{l + 6{y)r 



\Z — 7/1^^ 



that is exactly estimate 



6 Comments 

6.1 Well-posedness of the boundary layer system 

As mentioned several times in this paper, the well-posedness of system (ll.4p is not known 
without a structural assumption on uj, like periodicity or stationarity. We stress however that 
thanks to our estimates on the Green function G^^, the representation formula (j4.7|) defines a 
solution of (II. 4p for any C^'" boundary, cf the fourth section. Hence, the open issue is rather 
to find the appropriate functional space for uniqueness. 

Such difficulty does not arise when the Stokes operator is replaced by the Laplacian, or 
more generally by a scalar elliptic operator. Hence, one can show well-posedness in L°° of 

- Av = in 17^', v{y) = uj{yi) on d^^'' (6.1) 

if the function u) is bounded and Lipschitz. For the existence part, one may consider, for all 
n > 1, the solution v"^ of 

-At;" = in 17^' n Z)(0, n), v'\y) = uj{yi) on d (fl^^ D D{0, n)) . 

By the maximum principle, ||t;"||j;^oo < ||Li;||ioo, so that up to a subsequence it converges to 
some V in weak*. Straightforwardly, v satisfies (|6.ip . 
For the uniqueness part, let v G L°° satisfying 

-Av = in v{y) = on dO!'^ 

Let us show that w = 0. As we do not know the behavior of v at infinity, it does not follow 
directly from the maximum principle. In the case of a Lipschitz boundary w, we can conclude 
to the uniqueness in the following way. By the change of variables y := (yi, y2 — uj{yi)), the 
previous equation becomes 

div {A{y)Vv) = 0, 2/2 > 0, v\y^=o = 0, 

for some elliptic matrix A = {aij) with bounded coefficients. We extend A and v to {y2 < 0} 
by the formulas, f(yi,y2) := -f(yi,-y2), 

Om(yi,y2) := -am(yi, -y2), anj{yi,y2) ■= -anj{yi,-y2), 
aijiyi,y2) ■■= aij{yi, -y2) otherwise. 
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In this way, we get div (A{y)Vv) = on all M^. Harnack's inequality for elliptic equations 
(see [IH]) leads to 

sup {M + v) < C inf (M + v) 

\y\<R \y\<R 

for any R > and M such that M + u > 0. With M = max(0, — infu) and R going to 
infinity, we obtain that v = 0. 

6.2 Decay of correlations 

A key element in the paper is the estimate (14.40 on the Green function for the Stokes operator 
above an oscillating boundary. This estimate relies itself on the Holder regularity result (|5.4|) . 
In fact, with some more calculations in the same spirit, one could show a refined bound: 
for a C'^'" boundary uj and a source term / G C^'f^^D^ {0, 1)) (a,/i > 0), the solution of 
(I5.3p satisfies 

I|V-u'^||loo(£,£(o,i/2)) < (||'w||l2(£,£(o,i)) + ||/||cO.''(D'=(0,l))) ' 

with C independent of e. This yields an optimal bound on the Green function, that is 

E \dSGUz,y)\ + \VyPUz,y)\ < C^^^^. (6.2) 

\I3\<2 ^' 

Estimate (|4.1|) can in turn be improved as 

\v{uji,0,0) -v{uJ2,0,0)\ < Cn^^. 

Such bound is far easier to prove when p.4p is replaced by the scalar system (|6.ip . In this 
case, one may use a representation in terms of the standard two-dimensional Brownian motion 
B{m,t) = {Bi{m,t), B2{m,t)). If we denote (M,A^,/i) the probability space on which this 
Brownian motion is defined, 

v{io,0,0) = / {-io,0) {Bi{m,T{m))) dn 

where r is the exit time from Q^^uj) (see |25j). We now want to bound 

v{uji,0,0) - v{uJ2,0,0) = j (^uJi{Bi{m,Ti{m))) - uJ2{Bi{m,T2im)))^dfi 

where Tj is the exit time from Q^^uji). We remind that toi = lu2 over [— n, n], so the exit times 
are the same for brownian particles leaving in the region yi £ [—n,n]. Hence, 

|t>(u;i, 0, 0) — t;(a;2, 0, 0)1 < 2 max ||a;j||Loo /i (|i?i(-, rj)| > n) 

2 = 1,2 



< 2 ^fi {T^n < r_i) + fi{Tn< r_i; 

where we denote by T±n the first time for which Bi = ibn/2, and T_i the first time for which 
B2 = —1. It is well-known that the distributions of these hitting times are 

n f —n?\ , ,„ . ^ 1 /—I 
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A straightforward calculation provides 

which gives the result. 

6.3 Optimality of the decay rate 

Theorem H] shows that the boundary layer solution v converges at least as 1/2 ■ One may 
wonder if this result is optimal, that is if we can find roughness distributions for which the 

— 1/2 

speed of convergence is exactly given by 2/2 • Iii other words, is the constant Cj-q^o) of the 
theorem positive for some random distribution of roughness ? We have not so far been able 
to show optimality in this setting, but it can be established for the easier Dirichlet problem 

Au = 0, 2/2 > 
u = uj, y2 = 

where w is a given boundary data. Although simpler, this system shares many features with 
the original system p.4|) : 

• If w = uj{yi) is say L-periodic, the solution u{y) converges exponentially fast to the 
constant a := L^^ Jq uj{yi)dyi, as y2 goes to infinity. 

• liio belongs as before to the probability space (P, C, vr), one can show under assumption 
(H2) that 

y2^u{oJ,Q,y2) - a\'^ > ct^ > 0, a := E(a; a;(0)). 

Along the lines of [24^ pages 21-22], we will exhibit a stationary measure vr for which a > 0. 
Of course, vr is the law of the random process 99(0;, 7/1) := uj{yi), so that we just need to 
characterize the random initial data. Let G{uj,yi) a gaussian random process, of zero mean 
and covariance p{zi —yi), where p > is a smooth even function with compact support. Note 
that such process exists: take p = f * f, with / > an even smooth function with compact 
support. Then, its Fourier transform satisfies p = |/p > 0, which ensures the required 
positivity property 

c{zi)ciyi)p{zi-yi) = : [ | ^ c(zi)e^«^M'p(6 > 0. 

for any family c with compact support. Note moreover that this process defines almost surely 
smooth functions of yi: indeed, a simple calculation yields 

^^^\d'y^X{.,y,)fdyi^ = 2R{-lf p^^^\Q) < 00 

so that X{uj, •) is almost surely in the space Hj'^^{W) and therefore smooth. Finally, we 
introduce 

ip{uj,yi) = F{X{uj,yi)) 
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for a smooth increasing function F with values in (0, 1). We stress that ip satisfies (H2) as p 
has compact support. Wc will show that the corresponding a is positive. Suppose a contrario 
that fj = 0. For ?/2 > 1, wc introduce the measure tt^^ associated to the gaussian process with 
variance p{zi — yi) but mean 

m{zi,y2) = J p{zi-yi)g{yi,y2)dyi 

where g will be given later. Note that tt^^ is associated to the random initial data 

:= F{X((jj,yi) + 'm{yi,y2)). 

Standard computation yields 

R{i^,y2) ■■= = expN g{yi,y2)ip{uj,yi)dyi 

~\ j j P{^i-yi)9{zi,y2)g{yi,y2)dzidy^. 

and 

j \R{uJ,y2)\^d'K = e"^y^\ H{y2) '= j j P{zi - yi)g{zi,y2) g{yi,y2)dzidyi. 
If C7 = 0, then a simple Cauchy- Schwartz inequality 

j ^Jm\u{'^,^,y2)-a\d'K^^ < e5^(2'2) J y^\u{uj,0,y2) - al'^dir 

and goes to zero as 7/2 — ^ +00 if H is hounded from above. 

Let u and v solutions associated to the initial data ip{u!,yi) and (^^'^(a;, yi). As m > 0, by 
monotonicity, v > u. We can express u and v in terms of the Poisson Kernel, so that 

{v -u){uj,0,y2) = - 2^"! 2 ("^^"(^'^i) ~ 



Now, we define for 2/2 > 1 



9{yi,y2) = —=G[ — 

Vm \y2 



where G > has compact support, G = 1 over (—1, 1). On one hand, with this definition of 
g, one can check that 

sup H{y2) = sup ^ / / - yi)G ( — ) G (—] dzidyi < +00 

y2>i 2/2>i J J \y2j \y2j 



2/2 

On the other hand, one has 



/ Vmiv - u){uj,0,y2)dTr > C ( F'{X{uj,0))dTT] -^^^ ^m{yi,y2) dyi 

J Jr\J j yi + yi 



> C" I p{s)ds > 0. 



This implies that the quantity J (^(a;, 0, 2/2) — oi)diT does not go to zero, leading to a 
contradiction. 
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Appendix: Measurability of w 

We want to show here that 

/ Gu){z,y) eidy, z in , w{uj, z) := otherwise 



W[ljJ, z) 



defines a measurable function from P to Hl^J^). Let < ^ 1 a sequence of smooth 
functions with compact support, ^n\{-n,n) = 1- We define 



•Jn := / G^{z,y) {ifnei) dy, z in 0*'(a;), ^i;„(w,z):=0 otherwise. 

Jy2=0 

Note that Wn is the (unique) solution of 

' - Awn + Vg„ = 0, X G \ {y2 = 0}, 
div Wn = 0, xe 

_ [Wn\\y2=0 = 0, [d2Wn - {0,qn)]\y2=0 = (-<^n,0), 

satisfying J^tif^^^ |Vt(;„|^ < +00. By the dominated convergence theorem applied to the 
integral formula, we get that Wn ^ w in Lf^^. By the Cacciopoli inequality, the convergence 
is also true in Hj^^^. Thus, we just have to show measurability of Wn- 

Let us define 

V := {vG HHR^), div = 0} , := {v€ V, v\^2\^m^^^ = o} . 



29 



Following the lines of [HJ pages 15-16] it can be shown that the application uj vr(a;), where 
tt{u)) £ C(y,V) is the orthogonal projection from V to Kj, is measurable. Now, Wn is the 
unique fixed point of the contraction 

w 1-^ ^vr(tj) {w - Vn) 

where u„ is the unique function of H^{E?) satisfying 




The measurability of Wn follows. 
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